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Abstract 


Determining  optimum  distribution  of  fiber  through  rational  design  of 
corrugated  fiberboard  could  result  in  significant  reductions  in  fiber 
required  to  meet  end-use  conditions,  with  subsequent  reductions  in  price 
pressure  and  extension  of  the  softwood  timber  supply. 

A  theory  of  thin  plates  under  large  deformations  is  developed  that  is 
both  kinematically  and  physically  nonlinear.  This  theory  is  used  to 
analyze  the  strength  of  a  singlewall  fiberboard  element  which  is  modeled 
as  a  composite  plate  structure.  Either  an  elastic-type  instability  or  a 
compression  failure  mode  is  allowed.  The  results  of  this  analysis  are 
used  to  explain  how  to  adjust  facing  and  corrugating  medium  thicknesses 
to  maximize  the  edgewise  compression  strength  per  board  weight. 


Keywords:  Buckling,  compression,  corrugated,  design,  elasticity,  fiber- 
board,  model,  nonlinear,  optimize,  plates,  strength. 
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Introduction 
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In  the  design  of  fiberboard  containers,  one  of  the  basic  problems  is  to 
distribute  the  material  between  facings  and  corrugating  medium  to  maxi-  , 
mize  top-to-bottom  compression  strength.  Both  short-column  edgewise  .  1  -j  /  / 

compressive  strength  and  board  stiffness  contribute  to  the  overall  ^  J..  —  < 
strength,  and  have  been  investigated. 


One  of  the  earliest  theories  of  efficient  design,  prepared  by 
3/ 

Kellicutt  (4),-  was  to  maximize  strength  by  maximizing  the  board 
stiffness.  Kellicutt  advocated  locating  most  of  the  material  in  the 
facings  to  maximize  the  moment  of  inertia  and  the  bending  stiffness. 
He  also  developed  a  formula  for  calculating  the  moment  of  inertia  for 
a  corrugated  cross  section. 
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The  theory  developed  by  McKee,  Gander,  and  Wachuta  (7)  accounts  for  both 
short-column  strength  and  board  stiffness.  It  is  shown  that  the  optimal 
design  strikes  a  balance  between  these  two  considerations,  with  short- 
column  strength  playing' an  important  role  because  of  the  nonuniform  axial 
stress  distribution  across  the  fiberboard  panel.  Their  theory  is  based 
on  the  theory  for  the  buckling  of  thin  plates  modified  empirically  to 
match  experimental  data.  The  short-column  strength,  which  enters  the 
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3/  Underlined  numbers  in  parentheses  refer  to  literature  cited  at 
the  end  of  this  report. 
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formula  for  overall  strength,  is  obtained  by  direct  experimental 
measurement . 

Moody  (8)  provided  an  analytical  method  for  determining  the  short-column 
strength  based  on  local  buckling  of  the  fiberboard  elements.  His  theory 
agrees  reasonably  well  with  certain  experimental  data.  It  uses  results 
from  design  procedures  developed  for  airplane  structures  based  on  the 
fundamental  work  of  Cox  (1).  The  theory  allows  the  use  of  nonlinear 
stress-strain  data.  However,  the  method  of  calculation  does  not  lend 
itself  to  the  analysis  of  a  large  number  of  parameter  values. 

Zahn  (12)  has  reconsidered  the  foundations  of  the  local  buckling  theory 
and  provided  a  computer  program  to  analyze  a  large  number  of  parameter 
values.  This  analysis  is  restricted  to  linear  stress-strain  relation¬ 
ships.  Both  corrugating  medium  and  facings  must  have  the  same  elastic 
properties.  The  theories  of  both  Moody  and  Zahn  cannot  accomodate  an 
arbitrary  takeup  factor  in  the  corrugated  fiberboard  except  by  an 
arbitrary  approximation  of  the  geometry.  Zahn's  theory  allows  a  mode  of 
buckling  with  nodal  lines  diagonal  to  the  compressive  load.  This  mode 
is  not  included  in  Moody's  analysis  or  in  that  of  this  paper.  Some 
comparison  of  Zahn's  method  with  that  developed  here  is  made  in  the 
section  "Numerical  Results  and  Comparison  .  .  .  ." 

The  purpose  of  this  paper  is  twofold.  One  is  to  formulate  a  theory  for 
short-column  strength  of  the  local  buckling  type  which  allows  nonlinear 
elastic  behavior,  different  elastic  properties  of  medium  and  facing,  and 
an  arbitrary  takeup  factor,  and  which  provides  a  computer  program  which 
can  be  easily  used  to  generate  results  for  a  large  number  of  parameter 
values.  The  mathematical  analysis  is  summarized  and  an  iterative  method 
for  the  solution  of  the  buckling  equation  discussed.  An  empirical 
formula  is  used  to  fit  the  nonlinear  compressive  stress-strain  data,  and 
numerical  results  are  compared  with  experimental  data. 

Our  second  aim  is  to  examine  the  foundations  of  the  nonlinear  theory  of 
plates  on  which  our  buckling  theory  is  based.  A  theory  for  the  large 
deformation  of  plates  is  developed  which  allows  both  kinematical  and 
physical  nonlinearity  as  well  as  anisotropic  elastic  properties.  We 
know  of  no  other  theory  with  this  generality  which  is  needed  for  a 
complete  description  of  paper  behavior.  While  this  theory  is  not  used 
here  in  its  full  generality,  it  is  presented  as  a  basis  for  future 
contemplated  theoretical  and  experimental  research.  An  elastic  stability 
theory  is  developed  here  and  used  to  analyze  the  strength  of  a  short 
fiberboard  element  modeled  as  a  composite  plate  structure. 
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Short-Column  Compressive  St rength 


Buckling  Model  Development 

If  an  element  of  corrugated  fiberboard  is  under  a  uniform  compressive 
stress  applied  parallel  to  the  corrugations,  it  may  fail  due  to  local 
elastic  buckling.  This  mode  of  buckling  is  characterized  by  a  large 
number  of  lateral  deflection  half-waves  in  the  direction  of  the  load. 
Their  length  scale  is  that  of  the  microstructure.  We  note  that,  before 
local  buckling  can  take  place,  the  material  may  fail  due  to  the  ultimate 
compressive  stress  being  reached. 

To  develop  the  mathematical  model  for  local  buckling,  the  corrugated 
fiberboard  is  represented  as  a  composite  plate  structure  (fig.  1 1 .  The 
nonlinear  theory  of  plates  is  used  for  each  element  and  the  solution  of 
the  governing  differential  equations  for  local  buckling  is  discussed  in 
the  section  "Local  Instability  .  .  .  ."  The  most  serious  assumptions  in 
the  theory  are  that  the  plates  are  isotropic  and  that  linear  stress- 
strain  equations  can  be  used  with  an  additional  ad  hoc  assumption  which 
substitutes  the  tangent  modulus  for  Young's  modulus. 

A  buckling  equation  is  obtained  by  satisfying  boundary  conditions  that 
require  that  the  plate  elements  are  rigidly  fastened  and  that  the  junc¬ 
tures  be  in  moment  equilibrium.  The  buckling  equation  is  given  by 
(3.3.13)  in  the  form: 


K!  +  K2  +  K3  =  0 


(2.1.1) 


where  K.  (i  =  1,2,3)  are  plate  stiffness  coefficients  for  each  plate 
element  given  by  (3.3.16): 


K. 


E.t2  J  (l-V2)o. /E . 

_ l  l  v  l  i  l 

or.tan(a.b)  +  p.tanh(p.b) 
1  11  1  11 


(2.1.2) 


where  a.  and  p.  are  defined  by  (3.3.15).  The 
plate  element  (fig.  1),  cr  is  the  compressive 
tangent  modulus,  t.L  plate  element  thickness, 


subscript 
stress  at 
b.L  plate 


refers  to  .the 
buckling,  E^  the 
width,  and  L  the 
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Figure  1 . --Structural  model  for  singlewall  corrugated  fiberboard. 


length  of  the  sheet  in  the  direction  of  the  load.  In  eqs .  (3.3.15)  n  is 
the  number  of  half-wave  buckles. 

The  dimensional  parameters  in  fig.  1  can  be  derived  from  the  more  commonly 
understood  corrugated  fiberboard  geometry.  We  have 


hl  = 


h2  ♦ 


where 


h2  =  facing  thickness 

h^  =  corrugating  medium  thickness. 
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From  the  takeup  factor,  TF,  and  the  flutes  per  foot,  F,  we  obtain  the 
effective  spacing  between  flute  tips  as 


p  (in.)  =  1 2/F 


and  the  half-width  of  each  plate  element  as 


p2  (TF2- 1 )  -  4(hf  -  h3)2 
al=  8  p  (TF- 1 ) 

a2  =  p/2  -  a 

a^  =  p  x  TF/4  -  a^ 


where 


hj.  =  flute  height. 


It  is  found  that  the  available  stress-strain  data  for  paper  under  edge¬ 
wise  compression  can  be  fit  quite  accurately  up  to  the  ultimate  stress 
with  the  empirical  formula: 


a(e)  =  Cj  tanh(c2e)  +  c^e  (2.1.3) 

Here  O  and  e  are  axial  engineering  stress  and  strain  derived  from  the 
applied  load  per  undeformed  area  and  the  displacement  per  initial  length 
O  is  identical  with  N^/h  and  e  with  in  the  stress-strain  relation 

(3.1.25).  Details  on  how  coefficients  c.  are  determined  to  fit  the  data 

are  given  in  "Empirical  Stress-Strain  Formula."  A  curve  of  form  of 
eq.  (2.1.3)  is  fitted  to  data  for  facing  and  corrugating  medium  (fig.  2) 
The  tangent  modulus  is  given  by 


E(e) 


o'  =  CjC2/cosh  (c2e)  +  c3 


(2.1.4) 


Since  buckling  usually  occurs  when  the  strain  is  greater  than  the  propor 

tional  limit,  the  stress  a.  and  modulus  E.  in  (2.1.2)  are  taken  to  be 

l  i 


Figure  2. --An  empirically  derived  curve 

cj(e)  =  1,842  tanh  (196. 2e)  +  70,009  e  lb/in.2. 


given  by  eqs.  (2.1.3)  and  (2.1.4)  using  coefficient  values  appropriate 
to  either  facing  or  medium.  With  stress  and  modulus  given  in  terms  of 
strain  by  (2.1.3)  and  (2.1.4),  we  see  that  (2.1.1)  and  (2.1.2)  are  equa¬ 
tions  for  the  strain  £  at  buckling.  For  typical  fiberboard  properties, 
figure  3  shows  the  stiffness  coefficients  for  the  three  plate  elements 

(fig.  1)  and  their  sum  as  functions  of  strain  E.  The  number  of  half- 

4/ 

waves  is  taken  as  n  =  7.—  The  smallest  zero  of  the  sum  corresponds  to 
the  primary  mode  of  buckling  (fig.  4).  The  higher  zeros  correspond  to 
modes  in  which  the  plate  elements  undergo  wavier  deformations.  The 
primary  mode  that  has  the  smallest  strain  is  of  interest.  In  addition, 
n,  the  number  of  half-waves,  must  be  chosen  to  minimize  further  the 
primary  buckling  strain. 

The  analytical  solution  of  eq.  (2.1.1)  must  be  arrived  at  by  a  method 
of  successive  approximation.  It  is  found  that  the  following  algorithm 
converges  rapidly  to  an  accurate  buckling  strain.  If  the  buckling  is 
initiated  by  element  2,  start  by  solving  (2.1.1)  for  K^: 

k2  =  -«,<•  k3) 


4/  Additional  input  to  (2.1.1)  are  at  L  =  1.0,  t^  =  0.0194, 
t2  =  0.012,  t3  =  0.0074,  bj  =  0.0537,  b£  =  0.280,  b3  =  0.206,  and 

V1  =  V2  =  V3  =  °268- 
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Figure  3 . --Stiffness  coefficients  for  A-flute  construction  with  36  flutes  per 
foot,  0.18-inch  flute  height,  1.56  takeup  factor,  and  stress-strain  curves 
for  both  components  given  by  figure  2.  Facing  thickness  is  0.012  inch; 
medium  thickness  is  0.0074  inch. 
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Figure  4. --Primary  buckling  mode. 


Then,  solve  this  equation  for  tanCf^t^)- 


P2 

tan(a2t>2)  =  -  —  tanh(p2b2) 


E  t^ 
2  2 


a. 


\J (l-v^)o  /e 


(2.1.5) 


Next,  define 


Y  = 


/ 


12(l-vpa2/E2 

nnt. 


(2.1.6) 


Then,  (3.3.15)  becomes 


a2  =  n/i  Vy-1  ,  P2  =  nn  Vy+1 


(2.1.7) 


Introduce  (2.1.7)  in  (2.1.5)  and  invert  the  tangent  function.  After 
some  manipulation,  we  get 


Y  = 


+  - ^~~2  {" 

(nnb0)  '■ 


arctan[ 


tanh(nnb2VY+l) 


+ 


E2t2Y 

yi 2 ( Y- 1 )  (Kj+  K3)  j 


(2.1.8) 


Solve  (2.1.8)  by  iteration,  the  current  estimate  for  y  going  in  the  right 
side  to  generate  a  new  estimate  for  y.  The  steps  in  the  iteration  are 

1.  Initialize  n  as  n~L/b  where  L  is  the  length  and  b  the  width 
of  the  plate. 

2.  Initialize  y  =  yt  as  slightly  greater  than  1. 
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3. 


From  eqs.  (2.1.3),  (2.1.4),  and  (2.1.6), 


Cj  tanh(c2£)  +  c2£ 

2 

CjC2  sech  (c2e)  +  c3 


2  2  2 
n  n  t2 


12(1-v2) 


Y 


2 

1 


This  equation  is  solved  by  Newton's  method  to  give  an  estimate 
of  the  buckling  strain  £. 

4.  Compute  a^(£),  o3(e),  Ej(e),  E3(£)  from  (2.1.3)  and  (2.1.4). 
Compute  dj,  a^,  Pj,  p^  from  (3.3.15).  Compute  and 

with  (2.1.2). 

5.  Compute  right  side  of  (2.1.8)  for  new  estimate  y  =  y2- 

6.  If 


h 

h 


<  error  bound 


go  to  step  7.  Otherwise  go  to  step  3. 

7.  Increase  n  by  1 .  Go  to  step  3.  n  is  incremented  until  the 
minimum  buckling  strain  £  is  reached. 

To  speed  convergence,  an  overrelaxation  scheme  is  used  in  the  form 


new  y  =  v  Yj  +  0*°)  Y2 


where  0  <  a  <  1  is  a  convergence  factor. 

Both  eq.  (2.1.8)  and  the  corresponding  one  assuming  buckling  is  initiated 
in  element  3  are  used  in  the  program.  We  have  found  that  this  algorithm 
makes  possible  an  efficient  computer  solution  for  economically  reducing 
the  analysis  to  a  design  curve. 
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An  Empirical  Stress-Strain  Formula 


An  equation  of  the  form 


o(e)  =  Cj  tanhCc^e)  +  c^£ 

has  been  found  to  provide  a  good  fit  to  the  experimental  stress-strain 
data  digitized  from  a  load-deformation  curve  for  a  paperboard  specimen 
in  edgewise  compression.  If  there  are  n  data  points  (o. ,ei)  that  define 

an  experimentally  produced  curve,  the  coefficients  to  the  equation  may 
be  initialized  and  then  adjusted  to  improve  the  accuracy  of  fit. 

From  the  equation  for  a(e)  and  its  derivative 


O'  (£) 


c  c 
12 


cosh  (<  c) 


c  . 
j 


as  £  approaches  infinity,  the  limiting  value  for  o  (c)  becomes 


a 

n 


+  c„£ 

3  n 


which  can 
For  o' (e) 

limit  for 


be  used  to  estimate  c^  when  c^  is  known. 

the  limit  is  c^.  At  a  stress  near  the  ultimate  level  the 
o(£)  is  nearly  satisfied;  from  its  tangent  slope  calculate 


-  °n  ”  °n-l 
c3  "  e  -  e 


n 


'n-1 


Now  at  £  equal  to  zero 


o'  (0)  =  CjC2  +  c3 
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from  which  an  estimate  may  be  obtained  for 


c  c  =  — 
12  e. 


This  is  used  to  calculate  after  c^  and  c^  have  been  calculated  above. 

With  a  computer  program  the  coefficients  may  be  refined  to  provide  a 
more  closely  fit  equation.  Generate  the  array 


(1  -  e)c. 


(1  -  e)c. 


(1  -  e)c„ 


(1  +  e)c 


(1  +  e)c„ 


(1  +  e)c„ 


from  some  small  value  e  and  the  initialized  coefficients.  Select  one 
coefficient  from  each  column,  that  is  c  ,  c  ,  and  c  (27  possible  com- 

binations),  and  compare  the  fit  equation  to  the  original  data.  Define 
the  measure  of  fit  as  the  average  relative  difference  between  the  stress 
data  and  the  functional  evaluations  at  the  strain  data, 


n 


Nf . !  T  ... 

f  n  . _  i=l 


o. 

i 


o(e. ) 


o. 

i 


Then  let  the  coefficients  Cj,  c^,  and  c^  take  on  the  new  values  of  c^, 

c  ,  and  c  that  minimize  M,. 
v’  w  f 


A  load-deflection  curve  for  a  42-pound  linerboard  in  the  cross-machine 
direction  was  plotted  via  an  autograph  recorder  (fig.  2).  Eleven  points 
were  read  from  the  curve,  adjusted  to  appropriate  units,  and  used  to 
determine  the  coefficients.  The  fit  was  accepted  with  an  average  rela¬ 
tive  difference  of  1.02  percent.  This  curve  appears  as  number  5  in 
figure  5.  The  other  curves  in  Figure  5  are  determined  in  the  same  manner. 
Table  1  gives  the  coeficients  and  the  ultimate  stresses  for  these  curves. 
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1/  Average  coefficients  for  n  curves  are  found  to 
satisfy  the  following:  av.  C  =  (IC.)/n;  av.  C,  =  (XC 


Table  2. --Edgewise  compression  strength  analyses  compared  to  observations 


Flute  1 

size 

Flutes 

per 

Flute 

height 

Takeup 

factor 

Facing 

Thickness 

Measure- 

Stress-strain 
curve  number 
(fig.  5) 

Edgewise  compression 
strength 

Differ¬ 

ence 

Predicted 
f a i lure  *  ^ 

foot 

method 

Facing 

Medium 

Observed 

Predicted 

Type 

Component 

In . 

In. 

I*Ll 

Lb/ in . 

Lb/ in. 

Pet 

A 

(6) 

36 

0.  177 

1.55 

0.0118 

0.0073 

Stylus 

1 

13 

44.6 

B 

F 

(10) 

1 

2 

14 

13 

44.6 

46.7 

B 

B 

M 

M 

2 

14 

44.3 

B 

M 

45.8 

2/Ub.o 

-1.75 

A 

(6) 

36 

.174 

1.55 

0133 

.0075 

TAPP  I 

3 

15 

48.8 

B 

F 

3 

16 

46.5 

B 

M 

4 

15 

49.7 

B 

M 

4 

16 

46.5 

B 

M 

45.8 

47.9 

4.59 

A 

(6) 

36 

.  180 

1.55 

.0106 

.0080 

Stylus 

5 

17 

56.9 

B 

M 

5 

18 

52.3 

B 

M 

6 

17 

56.0 

B 

M 

6 

18 

51.5 

B 

M 

53.8 

54.1 

0.558 

A 

(5) 

36 

.  178 

1.55 

.0078 

.0084 

Stylus 

7 

19 

29.6 

B 

F 

7 

20 

29.4 

B 

F 

8 

19 

30.0 

B 

F 

8 

20 

26.8 

B 

F 

31.4 

29.2 

-7.00 

A 

(5) 

36 

.  178 

1.55 

.0185 

.0084 

Stylus 

9 

19 

59.8 

B 

M 

9 

20 

63.7 

B 

M 

10 

19 

63.5 

B 

M 

10 

20 

67.3 

B 

M 

65.5 

63.6 

-2.90 

A 

36 

.177 

1.53 

.0108 

.0080 

Stylus 

u 

21 

51.2 

B 

M 

11 

22 

49.0 

B 

F 

12 

21 

49.1 

B 

F 

12 

22 

68.1 

B 

F 

45.6 

49.3 

8.19 

B 

50 

.0961 

1.36 

.0108 

.0080 

Stylus 

11 

21 

55.3 

c 

M 

11 

22 

54.2 

C 

F 

12 

21 

56.0 

c 

M 

12 

22 

56.4 

c 

M 

48.5 

55.4 

14.3 

c 

42 

.  142 

1.46 

.0108 

.0080 

Stylus 

11 

21 

56.5 

c 

M 

11 

22 

55.3 

c 

F 

12 

21 

52.8 

B 

F 

12 

22 

51.7 

B 

F 

46.4 

54.1 

16.6 

1/  B,  buckling;  C,  compression;  F,  facing;  M,  medium. 
2/  Average  of  4  predictions. 
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Numerical  Results  and  Comparison 
with  Experimental  Data 

The  buckling  strain  solved  from  (2.1.1)  yields  the  stress  in  each  element 
from  the  respective  stress-strain  curve.  These  stresses  distributed 
over  the  combined  board  cross-sectional  area  define  the  failure  load  per 
unit  length  of  board  and  predict  the  specimen  edgewise  compressive 
strength . 

Comparisons  between  predicted  and  experimentally  determined  edgewise 
compression  strengths  are  presented  in  table  2.  Some  of  the  experimen¬ 
tal  values  are  taken  from  data  collected  for  other  FPL  investigations 
(5,6).  The  remaining  data  are  from  tests  performed  for  this  study.  The 
reported  observations  are  average  short-column  strengths.  All  available 
FPL  short-column  data  were  used  where  the  paper  stress-strain  data  were 
also  available.  The  predictions  are  derived  from  unpublished  load- 
deformation  curves  for  the  component  papers  of  the  short-column  specimens. 
F.q.  (2.1.3)  is  plotted  for  each  fit  to  the  paper  data  (fig.  5).  Two 
representative  curves  are  given  for  each  paper,  one  with  a  maximum  load 
above  and  one  below  the  average  maximum  load. 

All  materials  were  preconditioned  at  80°  F  and  30  percent  relative 
humidity  (RH)  and  conditioned  at  73°  F  and  50  percent  RH  prior  to 
testing  at  73°  F  and  50  percent  RH.  The  paper  component  were  tested 
for  edgewise  compression  in  the  cross-machine  direction  utilizing  the 
support  device  described  by  Jackson  (3).  The  rate  of  strain  was 
0.01504  inch/ inch/minute .  The  short-column  strength  values  were  then 
adjusted  to  the  components'  strain  rates  according  to  Moody  and 
Koning  (9).  The  construction  of  these  specimens  is  described  in 
TAPP  I  Standard  T811  (1J). 

The  predicted  location  of  failure  and  either  the  buckling  or  compression 
mode  are  identified.  Compression  failure  is  defined  as  occurring  when 
the  applied  load  produces  a  paper  stress  equal  to  a  component's  ultimate 
stress  before  buckling  can  occur.  One  observes  that  the  prediction  is 
sensitive  to  the  stress-strain  curve  shapes.  An  average  prediction  is 
thus  derived  from  the  permutations  among  curve  replicates. 

The  A  flute  data  shows  the  theory  to  agree  well  with  experiments.  A 
buckling  model  unadjusted  to  empirical  observations  would  normally  yield 
high  predictions.  It  is  difficult  to  design  a  model  that  accomodates 
all  possible  deformations.  It  appears  that  our  use  of  the  tangent 
modulus  theory  and  the  assumed  isotropic  behavior  compensates  for  this 
effect . 

For  the  B  and  C  flute  data  compression  failure  is  usually  predicted  at 
the  ultimate  stress  for  one  of  the  components.  The  compression  stress- 
strain  data  were  produced  with  the  paper  restrained  against  lateral 
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deflections.  The  accuracy  of  this  procedure  was  checked  against  tensile 
data.  The  stress-strain  curves  in  opposite  directions  were  found  to 
agree  (3).  It  is  logical  to  expect  a  plate  element  between  flute  tips, 
not  fully  restrained  against  lateral  deflections,  to  fail  in  compression 
at  a  stress  below  the  measured  ultimate.  The  B  and  C  flute  data  support 
this.  It  remains  to  be  investigated  how  to  best  adjust  the  theory  to 
utilize  the  stress-strain  relationship  for  predicting  compression  failure 
with  greater  accuracy. 

One  can  compare  our  analysis  to  Zahn's  by  altering  the  corrugated  fiber- 
board  geometry.  In  figure  1,  let  h^  be  dependent  on  2aj,  h^,  p,  and  TF. 

Then  let  2a^  vary  from  zero  to  some  small  length.  The  predicted  short- 

column  strength  varies  with  a  changing  geometry  (fig.  6).  The  point  at 
2a^  equal  to  0.0537  inch  is  compatible  with  nominal  A  flute.  Zahn's 

analysis  assumes  a  triangular  truss-core  shape  such  that  2a^  equals  zero. 

An  approximation  to  Zahn's  geometry  can  be  obtained  as  follows:  Let  2a^ 

equal  zero  and  reduce  2a^  by  2a  .  Calculate  the  buckling  stress  for 

this  new  structure.  From  this  stress,  determine  the  edgewise  compression 
strength  for  the  stress' distributed  over  the  original  geometry.  Zahn's 
method  yields  a  consistently  lower  strength  prediction  than  ours.  The 
reason  for  this  is  that  in  the  limit  as  a^  approaches  zero  our  mode  shape 

does  not  approach  the  correct  mode  shape  used  by  Zahn  for  the  triangular 
truss-core.  Zahn's  program  apparently  is  not  compatible  with  the  struc¬ 
ture  for  calculating  accurate  buckling  loads  for  fiberboard  structures 
with  takeup  factors  that  result  in  a  2a^  greater  than  zero. 

From  the  short-column  strength  prediction,  the  effect  of  distributing 
fiber  through  the  board  may  be  examined.  Figure  7  results  from  succes¬ 
sively  solving  eq.  (2.1.1)  over  an  array  of  facing  and  medium  thick¬ 
ness  combinations.  Solid  line  contours  of  equal  short-column  strength 
are  produced  via  an  interpolating  scheme  that  locates  all  the  facing  and 
medium  pairs  that  produce  equal  edgewise  compression  strengths.  Like¬ 
wise,  the  dashed  lines  define  the  board  weights  from  the  combined  paper 
weights  and  the  takeup  factor.  In  this  case  a  0.012-inch  paper  was 
assumed  to  weigh  42  pounds  per  thousand  square  feet  and  proportional 
equivalents  were  applied  to  other  thicknesses. 

As  there  is  yet  little  basis  for  selecting  a  stress-strain  curve  shape 
dependent  upon  paper  thickness  and  content,  average  stress-strain  curves 
were  determined  from  the  stress-strain  data.  The  average  linerboard 
curve  and  average  medium  curve  are  shown  (fig.  5);  their  coefficients 
are  given  in  table  1.  These  representative  curves  were  used  in  all  the 
calculations  to  represent  the  facing  and  medium  material  over  all  thick¬ 
nesses.  Additionally,  the  structure  analyzed  was  A  flute  with  the  same 
geometry  as  that  used  to  produce  figure  3. 
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Figure  7  identifies  regions  for  distributing  fiber  most  efficiently  by 
the  relative  slopes  between  strength  and  weight  contours.  As  an  example, 
given  the  board  with  a  0.018-inch  facing  and  a  0.00762-inch  medium,  the 
short-column  strength  is  read  at  65  pounds  per  inch  C 1 b/ in.).  For  a 
desired  strength  increase  of  5  lb/in.,  the  facing  thickness  may  be  main¬ 
tained  with  the  medium  thickness  increased  to  0.00817  or  the  medium  may 
be  held  constant  while  the  facing  is  increased  to  0.0198.  Both  alterna¬ 
tives  yield  a  70  lb/in.  short-column.  However,  the  combined  board  weight 
increases  from  168  to  171  pounds  by  adding  fiber  to  the  medium  and  to 
180  pounds  by  adding  fiber  to  the  facings.  It  is  thus  more  efficient  to 
add  material  to  the  medium  for  a  board  from  the  region  of  the  curve  with 
strength  contours  shallower  than  the  weight  contours. 

Where  the  strength  contours  are  steeper,  the  board  strength  is  most 
efficiently  improved  by  adding  fiber  to  the  facings.  This  becomes 
apparent  by  moving  from  a  133-pound  board  with  0.00966-inch  facings  and 
a  0.012-inch  medium  to  either  a  145-pound  board  with  the  same  facings 
and  a  0. 0142-inch  medium  or  a  139-pound  board  with  0.0104-inch  facings 
and  the  same  medium.  In  either  case,  the  short-coLumn  strength  increases 
from  50  to  55  lb/in. 

Where  the  strength  contours  become  tangent  to  the  weight  contours,  the 
optimum  design  is  achieved.  The  facings  and  medium  are  balanced  for  a 
desired  strength  with  minimum  fiber. 


Figure  7. --Edgewise  compression  strength  (lb/in.,  solid  lines)  and 
combined  board  weight  (lb/1,000  ft2,  dashed  lines)  over  an  array 
of  facing  and  corru  ing  medium  combinations.  The  stress-strain 
relationship  for  the  facing  and  medium  is  given  by  the  average 
curves  in  figure  5.  The  construction  is  A  flute  with  36  flutes 
per  foot,  0.18-inch  flute  height,  and  a  1.56  takeup  factor. 
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Equations  For  Large  Deflection  of  Anisotropic  Plates 


We  denote  Cartesian  material  and  spatial  coordinates  by  and 

The  locate  the  position  of  a  material  particle  in  the  reference 

state  while  locate  position  in  the  current  configuration.  Denote 

the  material  (or  Piola)  stress  tensor  (2)  by  P.  Component  P„  represents 

force  per  unit  reference  area  acting  in  the  j-direction  associated  with 
a  reference  face  perpendicular  to  the  i-direction. 


Because  tensor  P  is  associated  with  areas  and  directions  in  the  reference 
configuration,  it  is  a  convenient  measure  of  stress  to  use  for  the 
development  of  a  large  deformation  theory  of  plates.  For  no  body  forces, 
the  equations  of  equilibrium  can  be  written  as  (2) 

3P.  . 

=  0  (3.1.1) 


Note  that  x^  are  material  coordinates  and  that  eqs.  (3.1.1)  are  exact; 

they  contain  no  approximations  for  small  deformation.  Define  force  and 
moment  resultants  by: 
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where  the  integrations  are  taken  over  the  thickness  of  the  plate  and  h 
is  thickness  of  the  plate  in  the  reference  configuration. 

We  assume  that 


P31(h/2)  =  P31(-h/2)  =  0 

P32(h/2)  =  P32(-h/2)  =  0 

P33(h/2)  -  P33(-h/2)  =  q(x],x2) 


\  (3.1.3) 

J 


Integration  of  eqs.  (3.1.1)  and  use  of  conditions  (3.1.3)  yields: 
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>(3.1.4) 
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Next,  multiply  the  first  two  equations  of  (3.1.1)  hy  x.j  and  integrate. 
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Eqs .  (3.1.4)  are  the  equations  of  force  equilibrium  and  (3.1.3)  the  equa¬ 
tions  of  moment  equilibrium  of  a  plate  element. 

We  formulate  the  stress-strain  relations  through  use  of  the  work  princi¬ 
ple  : 


H  =  Pfc-F.. 

k  i  i  k 


(3.1.6) 


where  ff  is  the  strain  energy  density,  F.j  =  3f,./3x.  components  of  the 

deformation  gradient  tensor  and  the  dot  indicates  the  material  time 
derivative.  The  integrated  form  of  (3.1.6)  is: 


where 


J 


L. 


r  h/2 

I 

-h/2 

h/2 

-h/2 


P.  F.,  dx. 
k  i  i  k  3 


(3. 1.7) 


H  dx.j  is  the  plate  strain  enir«v  density. 
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To  proceed  further  we  make  the  usual  assumptions  for  a  thin  plate  con¬ 
cerning  the  form  of  the  displacement  components,  u.  =  t>.  -  x..  We  set 


u,  =  U  (x. ,x0)  -  x 


8w(Xl ,x2) 


1  “i  '■"i  >  2y  3  9x, 


8w(x  x  ) 

U2  =  U2(xl,x2)  ‘  X3  3^ 


U3  =  w(x1>x2) 


'N 


> (3.1.8) 


where  Uj,  U2>  and  w  are  components  of  middle  surface  displacement. 
The  components  of  finite  strain  are  given  by 


3u.  3u .  3u,  3u, 

2  E  =  _ -  +  _ 1  +  _ x  _ k 

ij  3x.  3x.  3x.  3x. 

J  1  1  J 


(3.1.9) 


Using  (3.1.8)  and  making  appropriate  thin  plate  approximations,  we  obtain 
from  Eq.  (3.1.9), 


E22  e2  "  K2X3 
E12  “  £12  ’  K12X3 


> (3.1.10) 


J 
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where  £ ^ ,  £^ ,  £^  are  middle  surface  strain  components  and  K ^ ,  K^,  kj2 
middle  surface  curvature  components. 
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>(3.1.11) 
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(3.1.12) 


We  assume  that  strain  energy  H  is  a  function  of  K ^ ,  F^,,  and  F.]:>. 


12' 


H  =  H(En,E22,E12) 


(3.1.13) 


With  (3.1.10),  we  see  that  the  plate  strain  energy  density  is  a  function 
of  £,  ,  £2,  c]2,  Kj ,  k2,  k]2. 

1  =  I(e1,£2,£]2;K1,K2,K12)  (3.1.14) 
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I 


With  form  (3.1.14)  and  (3.1.8),  the  work  energy  principle  (3.1.7)  yields: 


35  9U1  35  911 

Nn}  axj  +  (at^  "  N22)ax2 


1  35  3U1  1  3V  3U9 

2  a{.]2  N21}  3x2  2  9fcJ2  "  N12)ax] 


+  +  1  +  M  -  M  X  9V>’ 

la€j  3xj  2  ae12  8x2  N31  N13)ax] 


+  rJ*l  Jh?  +  1  Qw  -.aw 

vae2  3x2  2  ae]2  aXj  23  32^3x2 


,ai  a2w  ai  a2w 

+  (3k7  +  Mll}  r2  +  fe  +  M22)7T 

1  dXj  2  3x^ 


2. 

+  (dr~  +  +  M9i)  =  o 


12 


12  2i  axjax2 


-25- 


(3.1 . 15) 
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l 


to  obtain: 


Substitute  from  (3.1.19)  in  (3.1.4)^ 
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3  /»>  3w  ,  „  3w  \  .  3  /»,  3w  .  3w  ..  _  n 

axj  (N11  aw J  N12  3x2)  3x2  (N22  3x2  N21  3x^  q  "  (3  j  2Q) 


We  next  make  an  appropriate  plate  approximation  to  eq.  (3.1.13)  by  ex 
panding  in  a  Taylor  series  in  x  about  x„  =  0  and  keeping  the  leading 
terms.  Use  of  (3.1.10)  yields: 


H  -  H(e1,e2,e12)  +  9x^ 

0  "  0 


1  3  H 

X3  2  2 

J  z  axr 


(3.1.21) 


With  (3.1.21),  the  plate  strain  energy  density 


dx^  becomes 


T  _  .  „  ,  .  A  1  ,3  3  H 

X  h  H  (£j,£2,£j2)  +  24  ^  2 

8X3 


(3.1.22) 


Use  of  (3.1.10)  and  a  chain  rule  yields: 


32H 

3x2 


-  HnK1  +  H22k2  +  H33k12 


+  2H12KjK2  +  2Hj3KjK]2  +  2H23K2K12  (3.1.23) 
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where 


2  2 

1=0  H  =  9  H  t 

11  _  2  ’  12  Se.ae.  ’  ‘ 

d£  j  1  Z 


(3.1.24) 


With  (3.1.22)  and  (3.1.23),  (3.1.16)  and  (3.1.18)  yield 


.  3H_  _J1  3  "U  2  ^22  2 

h  BEj  24  h  ^  1  aej  2  ‘ 
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(3.1.25) 


M11  =  "  12  h  *H11K1  +  H12K2  +  H13K12^ 

M22  =  "  12  h  ^ H 1 2 K 1  +  H22K2  +  H23K12*  >(3.1.26) 

M12  +  M21  =  "  12  h  ^H13K1  +  H23K2  +  H33K12* 


The  basic  equations  are  (3.1.4)  ,  (3.1.4)2>  (3.1.11),  (3.1.12),  (3.1.20), 
(3.1.25),  and  (3.1.26).  These  are  fifteen  equations  in  fifteen  variables. 

For  an  isotropic  material,  H  is  a  function  of  the  invariants  of  the 
strain  tensor. 


H  =  H(I ,  II) ;  I  =  e1  +  e2>  11  =  6^2-6^ 


(3.1.27) 
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Linear  isotropic  elasticity  corresponds  to  the  functional  form: 


H  =  E  t2  E  tt 
2(  1-v2)  1+v 


(3.1.28) 


where  E  is  Young's  modulus  and  v  Poisson's  ratio.  From  eq.  (3.1.28) 
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>  (3.1.30) 


Eqs.  (3.1.25)  and  (3.1.26)  become 


Njl  =  ^  (£j+V£2) 
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N12  =  N21  ~  1+v  £12 


> (3.1.31) 
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(VKj+K2) 

12( 1-v2) 

X  M  -  - 

Eh3  „ 

21 


6(l+v)  12 
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For  a  linear  orthotropic  material  the  strain  energy  function  takes 
the  form 


E 

H  =  2(1-v'v2)  (El*VlE2> 


£ 

2(1-v^v2)  (e2+Vl£lel)  +  G  C12 


(3.1.33) 


where  v2Ej  =  VjE2. 
We  obtain 


3H 

3e 


1 


!-VlV2 


(C1+V2e2) 


M3. 1.34) 


l-VlV2 


<VlVE2> 


8H 


de 


=  2Ge 


12 


12 


J 


11  1-viv2 


‘22  1"VjV2 


V2^1 

H  —  u  —  4 

12  “  21  1_v1v2  ’ 


H13  =  H23  =  ° 


>(3.1.35) 

H,,  =  2G 
33  J 
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Eqs.  (3.1.25)  and  (3.1.26)  become 


hEj 

'll  =  l-V^  (C1+V2£2) 


hE 

N22  "  l-v1V2  ^viei+e2)>  N12  =  2Ghe12 


\ 


>(3.1.36) 


J 


Elh- 


M11  "  12(1-VjV2)  (VV2K2) 


E2h“ 


(V.K.+K.) 


*22  ~  12(l-V]\/2)  ''11  *2 


m  *  m  -  h^G 

M12  M21  ~  6  12 


>  (3.1.37) 


A  Small  Pert urb ation  Stability  Theory 

Assume  that  the  prebuckled  deformation  is  homogeneous  and  involves  no 
bending.  That  is,  ,  c^,  e12,  w*  Njj,  N22,  and  NJ2  are  constant.  The 

governing  equations  of  the  prebuckled  deformation  are: 


e,  = 


au 

dx 


11 


=  h 


1 

1 

9H 

a  e 


.  .s 

2  aXo 


,  au,  au9 
ei2  =  2  (a^  +  a^) 


i 


n  -  h  ?JL 

N22  h  3e„ 


>(3.2.1) 


N  -  m  _  I  h  _9H 

N12  “  N21  “  2  h  3e 


12/ 


These  variables  are  perturbed  by  additional  quantities  at  buckling  which 
are  assumed  to  be  small  and  the  equations  linearized  accordingly.  We 
shall  not  write  down  all  the  equations  governing  the  perturbation  but 
just  consider  what  happens  to  eq.  (3.1.20).  It  becomes,  with  q  =  0, 
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a2M 


n 


2 

a  m 


22  +  3  <M12tM21> 


3x, 


ax; 


axj3x2 


a2w 

^un* 2N 

°X1 


32w  +  n  3?w  =  o 

12  ax:ax2  22  3x2 


(3.2.2) 


Use  (3.1.12)  and  (3.1.26)  to  write  the  moments  in  (3.2.2)  in  terms  of 
buckled  displacement  w. 


h_ 

12 


~4  a4  <*4 

Hn  +  H?2  4  +  (2H  +H  ) 

11  axj  22  ax4  12  33  ax2ax2 


+  2H 


a4w 


13  3„ 

3xj3x2 


+  2H 


fl4 

a  w 


=  N..  —  +  2N 

11  axj 


23  ax^x3 


a2w  +  N  a^w 

12  9xi9x2  22  ax2 


(3.2.3) 


In  this  equation,  w  is  the  displacement  due  to  buckling,  while  N.^,  N.^, 
and  N22  are  associated  with  the  prebuckled  deformation.  are  functions 

of  e  ,  e2,  and  eJ2  associated  with  the  prebuckled  deformation. 

For  the  linear  isotropic  material,  the  H. .  are  given  by  (3.1.30)  and 
(3.2.3)  becomes  1J 


Eh3  „2„2  _  „  a2w  .  „„  32w  .  „ 

- —  V  V  w  =  N  — r  +  2N.»  a"~  +  N 

i2(i-v2)  11  axj  12  9xi8x2 


2 

aw 

22  _  2 
3x2 


(3.2.4) 


where  V2w  =  32w/3x2  +  a2w/3x2. 
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;  i 


For  the  linear  orthotropic  material,  the  H..  are  given  by  (3.1.35)  and 
eq.  (3.2.3)  becomes 


D,  +  2(D'+C)  — +  D,  ^ 
9Xj  3x 3xJ 


a2  2  .2 

n  °  +  om  3  w  3  w 

n  ax2  12  3xi9x2  *22  ax2 


(3.2.5) 


where 


12(1-v1v2) 


V2°l  =  V1D2 


2  "  12(1-v1v2) 


(3.2.6) 


Note  (3.1.37)  can  be  written 


Mn  =  -  DjOCj+VjKj) 

M22  =  ‘  D2(V1K1+K2) 
M12  +  M21  =  "  2CK12 


(3.2.7) 


Local  Instability  of  Corrugated  Fiberboard 

The  preceding  stability  theory  is  applied  to  the  buckling  of  corrugated 
fiberboard  subject  to  axial  compression.  As  modeled  (fig.  1),  it  is 
a  periodic  structure  of  three  rigidly  connected  plate  elements,  with 
thicknesses  h^  and  widths  2a^.  Local  coordinates  are  chosen  for  each 

plate  element  with  the  origin  at  the  center  (fig.  8).  Let  the  prebuckled 
loading  be  an  axial  strain  e  in  each  plate  in  the  y-direction.  If  the 
Poisson's  ratios  of  the  material  of  each  plate  are  equal,  their  strains 
in  the  x-direction  will  be  the  same.  Hence,  the  plates  will  deform 
without  bending  and  form  a  figure  which  is  geometrically  similar  to  the 
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Figure  8. --Local  coordinates  for  plate  elements  of  model. 
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undeformed  figure.  We  assume  that  the  prebuckled  deformation  has  this 
feature  even  though  the  deformations  are  large  enough  to  necessitate  use 
of  a  tangent  modulus  theory.  Under  this  assumption  the  prebuckled 
deformation  is  homogeneous  so  that  the  equations  of  the  preceding  section 
can  be  used  for  each  plate  element. 

We  now  make  two  assumptions.  We  assume  that  each  element  of  the  model 
can  be  treated  as  an  isotropic  plate  and  use  eq.  (3.2.4).  Because  it 
is  found  that  buckling  takes  place  at  strains  beyond  the  linear  elastic 
limit,  we  account  for  this  by  using  the  tangent  modulus  for  E.  E  is 
therefore  a  function  of  the  strain  e.  It  is  the  cross-machine  modulus 
of  the  paper.  Paper,  of  course,  is  highly  anisotropic  so  that  a  more 
accurate  theory  would  use  (3.2.5)  rather  than  (3.2.4).  The  magnitude  of 
the  error  in  this  approximation  remains  to  be  investigated.  Also,  the 
tangent  modulus  assumption  which  is  ad  hoc  should  be  justified  by  showing 
that  there  is  a  strain  energy  density  H  which  leads  to  eqs .  (3.1.30) 
with  constant  E  replaced  by  the  tangent  modulus  E(e). 

For  axial  compression,  N12  =  =  0,  N22  =  N  and  eq.  (3.2.4)  becomes 

V2V2w  +  12(1~V  )N  ^  =  0  (3.3.1) 

EhJ  By 


The  buckling  load  N  is  the  smallest  eigenvalue  for  which  eq.  (3.3.1)  has 
a  non-trivial  solution.  The  buckling  modes  are  assumed  to  take  the 
following  separated  form: 


w  =  f(x)sin(\  y)  ,  \  =  r— 

n  ’  n  L 


(3.3.2) 


which  corresponds  to  taking  the  edges  of  the  plate  elements  simply 
supported  at  y  =  0,  L.  Substituting  (3.3.2)  in  (3.3.1)  yields  the 
following  equation  for  f(x): 


2  2 
ft  f 


**  f  =  0 

Eh3  " 


(3.3.3) 
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The  general  solution  of  (3.3.3)  is 


f (x)  =  CjCoshto^x)  +  c2cosh(a2x) 
+  c^sinhCofjX)  +  c^sinh(c<2x) 


(3.3.4) 


where 


“l 


+  slng^S 

"V  A2Eh- 


)N 


(3.3.5) 


The  primary  buckling  mode  is  even  in  x  in  each  plate  element  (fig.  4). 

Hence,  c„  =  c,  =  0  in  solution  (3.3.4)  and  c,  and  c„  are  to  be  determined 
3  4  1  2 

from  appropriate  boundary  conditions.  The  solution  in  each  plate  element 

is 


f.(x)  =  c.,cosh(a..x)  +  c . Ocosh(a . „x) 
i  ll  ll  i2  i2 


(3.3.6) 


where  index  i  denotes  the  element  and 


a., 

i  lX} 

“i2 


(3.3.7) 


For  boundary  conditions  we  require  zero  deflections  at  the  junction  of 
plate  elements, 


w j  (a  j )  =  w2^"a2^  =  W3^‘a3^  =  0  (3.3.8) 

equal  rotation  of  each  element  at  junctions, 

3w  (a  )  9w  (-a  )  3w  (-a  ) 

-ir - V2-  =  -V-  =  Rsin<V>  <3-3-9> 
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and  moment  equilibrium  at  junctions, 


-  H2x<-2>  -  H3x<-3>  =  0 


(3.3.10) 


Conditions  (3.3.8)  are  approximations  to  the  more  accurate  conditions 
that  would  require  equal  deflection  of  elements  at  junctions.  Conditions 
(3.3.8)  result  from  assuming  that  in-plane  buckling  displacements  in 
each  plate  element  are  negligible  compared  to  transverse  displacements 
.  Because  of  this  assumption  it  is  not  possible  to  impose  force 

equilibrium  at  each  joint.  This  assumption  has  been  made  and  justified 
for  aircraft  structures  (1). 

Eqs.  (3.3.8)  yield 


cosh(or  2^) 

il  cosh(a^a^) 


(3.3.11) 


while  (3.3.9)  and  (3.3.11)  combine  to  give 


Rcosh(of  . a  ) 

^  _ _ _ _  / ^  121 

12  a12cosh(a11a1)sinh(a12a1)  -  ^cosh^^a^  )sinh(a1  aj)  ^  * 


with  similar  expressions  for  c  „  and  c„  .  Using  (3.1.32),  (3.3.11),  and 
(3.3.12)  yields:  Z  J 


Ejh?  (°'ll‘°'l2^  R  sin(V) 

12(l-vp  '11tanh(°'11a1)  -  a12tanh(a12ai) 


with  similar  expressions  for  M2x(-a2)  an<^  ^3x^"a3^'  ^he  moment  equi¬ 
librium  equation  (3.3.10)  yields  the  buckling  equation: 


K1  +  K2  +  K3  =  ° 


(3.3.13) 
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where  the  K.  (i  =  1,2,3)  are  plate  element  stiffness  coefficients. 


K.  = 


Eihi  A/EiV.-v2) 


1  -  L'5a.2tanh(a.2ai)  +  L'V1tanh(a.1a . ) 


(3.3. 14) 


Write 


o.  =  N./h.  ,  t.  =  h./L  ,  b.  =  a./L 

liiii  ii 


a  =J- 1  La .  „  =  nn 

i  i2 


y^d-v^a./E. 


*\ 


n/tt . 

l 


-  1 


>  (3.3.15) 


Pi  =  Laii  =  nn\J 


yi2(l-Vj)0-/E. 


nnt . 
i 


+  1 


J 


Then  the  stiffness  coefficients  can  be  written  in  the  form 


E.t^  yo-vfo./Ej 

i  a.tan(a.b)  +  p.tanh(p.b.) 
i  ii  l  i  l 


K.  = 


(3.3.16) 


Because  a.  and  E.  are  functions  of  the  axial  strain  £,  eq.  (3.3.13)  is 

l  l 

an  equation  for  the  buckling  strain  £  on  the  structure.  Its  solution 
is  considered  in  the  section  "Buckling  Model  Development." 
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